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ABSTRACT 

We present a catalog of near-infrared (NIR) spectra and associated measurements for 886 nearby 
M dwarfs. The spectra were obtained with the NASA-Infrared Telescope Facility SpeX Spectrograph 
during a two-year observing campaign; they have high signal-to-noise ratios (SNR > 100 — 150), 
span 0.8-2.4 /rm and have R ^ 2000. Our catalog of measured values contains useful Tgff and 
composition-sensitive features, empirical stellar parameter measurements, and kinematic, photometric, 
and astrometric properties compiled from the literature. We focus on measures of M dwarf abundances 
([Fe/H] and [M/H]), capitalizing on the precision of recently published empirical NIR spectroscopic 
calibrations. We explore systematic differences between different abundance calibrations, and to 
other similar M dwarf catalogs. We confirm that the M dwarf abundances we measure show the 
expected inverse dependence with kinematic, activity, and color-based age indicators. Finally, we 
provide updated [Fe/H] and [M/H] for 16 M dwarf planet hosts. This catalog represents the largest 
published compilation of NIR spectra and associated parameters for M dwarfs. It provides a rich and 
uniform resource for the nearby M dwarfs, and will be especially valuable for measuring Habitable 
Zone locations and comparative abundances of the M dwarf planet hosts that will be uncovered by 
upcoming exoplanet surveys. 

Subject headings: stars: low-mass—stars: fundamental parameters—techniques: spectroscopic—stars: 
abundances 


1. INTRODUCTION AND BACKGROUND 

Recent years have seen a crescendo of observational 
and theoretical interest in the population of nearby M 
dwarfs. These stars, which make up t he majority of th e 
stars in the solar neighborhood (e.g. iReid et ^l2002h . 
have historically been difficult to observe due to their 
optical faintness. During the past two decades, great 
strides have been made in the refinement of techniques 
for identifying nearby M dwarfs. Multip l e groups have 
used parallaxes (e.g 


(2011 I.Tao et all mil 


M dwaris. Multip l e groups have 
Henry et al.l l2006t [Winters et al.l 


iRiedel et al.l 1201^ photome- 
try, and proper motions to ide ntify thousands of M 
dwarfs in the Solar neighborhood (IlI pine fc S harall2005l: 
Lmine fc GaidosI 120111: iLfeine et al.l l201,8l iFrith et al.l 


201,81: iGaidos et al.l l2014fl and througho ut the Galaxy 
(IWest et al.ll201lHDeshDande et al.ll2013fl . 

In parallel, tools for observationall y characterizing 
these stars have improved dramatically (iReid et al.l 
Hawley et alJ llOOd iGizis et HI 120021 : IReid et alT 
Leoine et al.l l200ffl . Even so, the precise measure¬ 


ment of their stellar parameters remains challenging. 
Precise stellar radii and masses have only been mea¬ 
sured fo r a small set pf ta r gets through interfer ome- 
try le.g. iDemorv et aH 120091 : iBovaiian et al.l 12012) and 
the study of eclipsing binaries le.g. iMorales et al.ll2009l : 
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iCarter et al.ll20ll]l . The measurement of M dwarf stel¬ 
lar compositions has also presented a formidable chal¬ 
lenge, as the mol ecular species in th eir atmospheres mask 
the continuum ([Gustatssonl 11989( 1. typically rendering 
the direct spectroscopic modeling techniques intractable 
(altho ugh progress is being made, e.g. lOnehag et al.l 
(|2012(l l. Instead of direct modeling, the detailed analysis 
of atomic line and molecular absorp tion in M dwarf spec¬ 
tra has also been deinonstr ated (jWoolf fc WallersteinI 
(200^ (200^ [Woolf et al.ll20nll . 


Rather than depending solely on the M dwarf spec¬ 
tra to derive abundance information, several groups have 
instead relied on abundance measurements of associ¬ 
ated stars. These calibrator systems typically contain 
an FGK dwarf in a common-proper-motion system with 
an M dwarf, separated by several arcseconds or more. 
The composition of the higher-mass star is known or 
readily determined using established spectr oscop i c tech - 
niques (as deployed for M subdwarfs in iGizi 
The M dwarf is then assumed to have the same com¬ 
position, and their photometric and spectroscopic prop¬ 
erties are used to construct empirical models that can 
then be applied to single stars. This empirical calibra¬ 
tion strategy was first implemente d for optical photom¬ 
etry iy — Kg, iBonfils et al.l 1200,^ . taking advantage of 
the strong and weak [Fe/H] dependenc e for M dwarf 
y arid Ar..i-band absorption, respectively (iGhabrier et al.l 
120001: iDelfosse et al.ll2Mot IBonfils et al.ll2005ll . Multiple 
groups have developed similar photom etric calibrations 


I 

2012 

1 

2014 


Heiazi et al .1 1201,5(1 . and iNeves et ahl (|2018aL I2(H^ have 


further leveraged these photometric [Fe/H] calibrations 
into more precise optical spectroscopic calibrations. 
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This empirical calibration technique has also been im¬ 
plemented for near-infrared (NIR) spectra, taking ad¬ 
vantage of the relative brightness of M dwarfs in the 
NIR. The fi rst of these was the RTs-ba nd calibration de¬ 
veloped by IRoias-Avala et afl (|201Clll . which leveraged 
the A"s-band Na i doublet and Ca i triplet, along 
with H 2 O absorption indices to accoun t for Tpff ef¬ 
fects. This technique was e xpanded bv iTerrien et afl 
(l2012bL , T12), r efined in iRoias-Avala et al.l (|2012ll 
and iNewton et a hl (120141) (h ereafte r N14), and further 
generalized in iMann et 51 (j2013a[l (hereafter M13a). 
Presently, the NIR calibration with the largest number 
and diversity of ca librators is that of M13a for early- 
mid M dwarfs, and IMann et al.l (I2014H (hereafter M14) 
for late-type M dwarfs; these calibrations provide [Fe/H] 
and [M/H], are well-calibrated across a wide range of 
[Fe/H] (—.8 < [Fe/H] < .5) and spectral type (late K to 
M8), and have good precision (^ 0.10 dex in [Fe/H] and 
[M/H]). 

In addition to the empirical techniques for mea¬ 
suring stellar composition, various groups have devel¬ 
oped similar techniques for the measurement of spec¬ 
tral type, effective temperature, stellar mass, and stel¬ 
lar radius, based on spectral type standards, stars 
with interferometrically-determined parameters, or stars 
with parameters measured through other means. No¬ 
tably, techniques applicable to high resolution opti¬ 
cal spectra froin HAR PS (iNeves et al.l l2013al [2014 
iMaldonado et all 1201511 and Keck ( Pineda et al.l l2013fl 
have achieved good precision for the M dwarfs observed 
in these programs. Techniques based on lower-resolution 
(R ^ 2000) NIR s pectra, using Teff-sensitive indices 
IjMann et al.l 12013^ (hereafter M13b) and atomic ab¬ 
sorption features ( Newton et al.l I2015H have been also 
been developed. The use of broadband photometry 
has also enabled precise measure ments of the tem - 
peratures and radii of these stars (IMann et al.l l2015ll . 
The resulting suite of techniques has been used ex¬ 
tensively in t he characterizat i on of Ke pler low-mass 
plane t hosts (iM uirhead et al.1 l2012alibl: l.lohnson et al 


2?iTa IHallard k all 1201 31: IMuirhead et all 1201 31 l~ 


Newton et al.ll201^ . 


The tools and data described above, combined with 
advances in structural (Baraffe et al.l 119981 : iDotter et al 


2007 ; iSoada et al.l 12013 ) and atmospheric (|Allard et al. 


201211 models, provide a fertile ground for a vari¬ 


ety of topics relat ed to M dwarfs , incl u ding Galac- 
tic evolution (e.g. 
lar s t ructure (e.g. 


Bochanski et alJ l2 01(1 Iml. stel- 

Lopez-Moralesl l2007t iReiners et ^ 


l2012t iFeiden fc Chaboved l2012t lTorresH2013ll . and ex- 
oplanets. Nearby M dwarfs are appealing targets for 
Doppler Radial Velocity (RV) and transit exoplanet 
searches because their low masses and luminosit ies make 
plane t s orbiting in their Habi t able Z ones (HZ iKastind 
Il993t iKopparapu et al.1 I2013L 1201411 more easily de¬ 
tectable than those around higher-mass stars. More¬ 
over, given the prevalence of M dwarfs, their plan¬ 
ets may represent the most common exoplanetary en¬ 
vironment in the Galaxy. Upcoming NIR ground-based 
RV surveys such as the Habitable-zo ne Planet Finder 
on th e Hobby-Eberly Telescope CHPF. iMahadevan et al.1 
IMll . CARMENES (iQuirrenbach et al.11201211 a t Caly 
Alto, SPIrou on the CFHT ( Thibault et al.l 1201211 . 
IRD on Subaru ([Kotani et al.l 1201411 . and iLocater 


on LBT as well as photome tric surveys like TESS 
([Ricker et al.1 l2015ll. PLATO (iRauer et al.11201411 . and 
MEarth ( Nutzman fc Charbonneaul 1200811 will specifi¬ 
cally target planets around nearby M dwarfs. These sur¬ 
veys will provide unprecedented sensitivity to planets or¬ 
biting the nearby M dwarfs, whose optical faintness has 
made them difficult targets for many existing RV and 
transit surveys. Despite being optimize d for optically 
brigh t FGK stars, results from Kepler ([Borucki et al.1 
l20inH clearly indicate that M d warfs host a large popula¬ 
tion o f small planets in the HZ ([Dressing fc Charbonneail 
I 2 OI 5 II . The characteristics of planets around M dwarfs 
will be useful for test ing theories of planet formation (e.g. 
iLaughlin et al.ll2004ll in the regime of low stellar and pro¬ 
toplanetary disk masses. 

In order to construct a well-characterized target list for 
HPF, we have carried out an extensive spectroscopic sur¬ 
vey of 886 nearby M dwarfs, using the SpeX spectrograph 
at the NASA Infrared Telescope Facility (IRTF). We take 
advantage of the NIR spectroscopic techniques described 
above to measure the metallicities of these stars, in order 
to form the largest catalog yet of nearby M dwarf compo¬ 
sitions based on spectroscopy. We also provide measure¬ 
ments of the Teff-sensitive indices, which can be used to 
estimate stellar mass, radius, or luminosity. This work 
echoes the spirit of previous large-scale effort s to charac¬ 
terize our most common stellar neighbors le.g.lReid et al.1 
119951 : iRlaz et al.ll2006l : iDeshoande et al.ll20i3l N14b 


In Section [2] we describe our observations and 
data reduction, including target selection and cross¬ 
identification with other catalogs. In Section [3| we de¬ 
scribe our spectral measurements and error estimates. 
In Section [5| we provide a brief summary and discussion 
of our measurements and comparisons to other catalogs. 


1.1. Published Results from this Sample 

This dataset has already enabled a variety of studies 
of the nearby M dwarfs. T12 developed a new H-band 
[Fe/H] calibr ation for M dwarfs, e xtending established 
techniques (iBonfils et al.1 120051 : iWoolf fc Wallersteliil 
I 2 OO 6 I : IRoias-Avala et al.l 1201011 and providing a valuable 
comparison for the development of improved calibrations 
(M13a) and a starting point for work with higher res- 
olution spectra (e.g. SDSS-HI APOGE E ! or IGRINS, 
[ Peshoande et al.1120131: iPark et al.ll2014D . ITerrien et all 
( 2012all employed the T12 [Fe/H] measurement tech¬ 
nique to characteri ze the low-mass eclip sing binary GM 
Draconis (CM Dra iMorales et all l2009ll . helping to re¬ 
move the [Fe/H]-related degeneracy and to clarify the 
disagree ment between models and observations of this 
system. ITerrien et al.l ([201411 again applied the T12 
[Fe/H] measurement technique to confirm the first low- 
mass members of the near by Coma Ber e nices cluster 
([Trumpled 1193811 . Finally, ITerrien et al.1 (1201511 lever¬ 
aged M dwarf activity catalogs (e. g. iGizis et al] 120021 : 
iRiaz et ahl 120061 : iGaidos et al.ll20l3 l and improvements 
in parallax measurem ents of the nearby M dwarfs (e.g. 
iDittmann et al.ll20l4 and the Research Consortium on 
Nearby Stars, RECON^) to develop new spectroscopic 
indices for M dwarf luminosity, surface gravity, and pos¬ 
sibly a-enrichment. 

® http://www.recons.org 





















































































































































A Near-Infrared Spectroscopic Survey of 886 Nearby M Dwarfs 


3 


2. OBSERVATIONS AND DATA REDUCTION 

Over a period of two years, beginning in 2011 May and 
concluding in 2013 August, we observed 886 M dwarfs 
with the SpeX NIR spectrograph on the NASA-IRTF 
(|Ravner et alJl200^ . This set of observations includes 
21 stars observed multiple times. 

2.1. Target Selection 

We drew M dwa rf targets pri marily from the 
iLepine &: Sharal (|2005L LSPM-N) and iLeoine fc GaidosI 
( 20 111 LGlll proper-motion based catalogs of nearby M 
dwarfs, and from published lists of M dwarf planet host 
and k nown wide binaries (drawn from iRoias-Avala et ahl 
[Moll that serve as abundance calibration stars. The 
majority of target stars (829) were taken from LGll, 
the remainder being the aforementioned binary calibra¬ 
tor stars, planet hosts, or that were t argeted in a con¬ 
temp oraneous high-resolution survey ([Deshoande et all 
1201311 . For stars from LGll and LSPM-N, we prioritized 
targets with redder optical-to-infrared color (V — J > 4), 
a good indication of mid-late M spectral subtype (M4 
or later). The resulting compilation of targets included 
both high-priority targets (calibrator stars, planet hosts, 
late-type M dwarfs) and low-priority targets (early type 
M dwarfs). 

From this super-set of potential targets we constructed 
our nightly observing plans, and during each night we ac¬ 
tively added and removed targets to account for changes 
in observing conditions. Our primary goal was to en¬ 
sure very high signal-to-noise (SNR > 200 per pixel) for 
high-priority targets, and our secondary goal was to ob¬ 
serve as many M dwarfs as possible with reasonable SNR 
(> 100). The typical resulting observing pattern was a 
series of groups (in position and time) which contained 
a few faint late-type (M4 or later) M dwarfs and several 
brighter early (earlier than M4) M dwarfs. Each group 
was sufficiently localized that a single nearby standard 
star (within one hour of time and 0.1 sec z of target star) 
could be used for telluric correction, and the same arc 
lamp exposure could be used for wavelength calibration. 
This strategy enabled a large number of M dwarfs to 
be observed, with minimal time losses from large slews, 
wavelength calibrations, and standard star observations. 

2.2. Data Reduction 

We obtained uniform data for all observations, op¬ 
erating SpeX in the short cross-dispersed (SXD) mode 
with the 0.3 x 5" slit, which produces R ^ 2000 spectra 
from 0.8 - 2.4 /xm. As described in T12, we extracted 
thes e spectra with the fac ility-provided SpeXTool pack¬ 
age (iCushing et al.l 1200411 . We tellmic-co rrected using 
the xtellcor program (jVacca et al.l 120031) and observa¬ 
tions of an AOV star (or similar), at an airmass and ob¬ 
servation time within 0.1 (sec z) and 1 hour of the target, 
respectively. The xtellcor program also flux-calibrates 
the spectrum, using the optical magnitudes (from SIM- 
BAD) of the stand ard star to es t imate its flux as a func¬ 
tion of wavelength. iVacca et al.l (j2003f l estimate that the 
flux calibration using this method agrees to within a few 
percent of fluxes derived using broadband magnitudes. 

® e.g. http://exoplaiiets.org 


A typical example spectrum is shown in Figure [T] All 
observations presented here took place prior to the SpeX 
detector upgrade, which occurred in the spring of 2014. 


2.3. Cross-Identification with other catalogs 


We cross-identified our targets with several published 
catalogs, in order to provide supplementary data for our 
target stars and to facil itate the use of our ca talog, Each 
target has a 2MASS ([Skrutskie et al.l I2006D identifier, 
and we use the 2MASS J2000 coordinates in our final 
catalog. The 2MASS ID also provides the most conve¬ 
nient uniform identifier for all our targets. We include 
in our catalog the 2MASS J, H, Kg photometry and the 
respective quality flags for each star. 

We cross-identified our targets with the UCAC4 
(jZacharias et al.ll2013[l catalog as well. UGAG4 provides 
optical magnitudes based on the AAVSO Photometric 
All-Sky Survey (APASS). Optical magnitudes are useful 
for M dwarfs as their optical-to-infrared colors are ex¬ 
cellent indicators of spectral type and Teg; we include 
UCAC4 aperture (A)-magnitucles (for 874 stars) and 14- 
magnitudes (for 732 stars) as indications of their optical 
brightness. 

Eor the majority of our targets, we include proper 
motions as measured by LGll or LSPM-N. For tar¬ 
gets without proper motions in these catalogs, we in- 
clud e measurements fro m UGAG4 or the PPMXL cata¬ 
log (jRoeser et al.ll201fll l. 

We also provide parallax measurements for a nearly 
half of our targets (427). Many of these parallax 
measurements are drawn from the REGONS database 


llao et al.l 120051 : iRiedel et al.l 12014 120101: iHenrv et al.l 
120061) and the newly-publis hed MEarth parallax cat¬ 
alog (jDittmann et al.l 12014 ). We also found par- 


19951 Ivan Leeuwenl 120071 IGould & Chaname 

2(1(11: 

Anglada-Escude et al.l 120121 Harrington et al. 

1993: 

Gatewoodll2008l IMonet et al.lll992l ILeoine et all 

200^ 


3. SPECTRAL MEASUREMENTS 

We describe here the various spectral measurements 
we performed on the M dwarf spectra. 


3.1. Band Indices and Spectral Type 

The NIR spectra of M dwarfs are dominated by the 
broad absorption of H 2 O, the strength of which is known 
to correlate strongl y with Teff and spectral type (e.g. 
iWilking et aLllTO^ . Using R ^ 2000 NIR spectra like 
ours, several groups have developed H 20 -based and other 
indic es that are well-co r related with spectral type and 
Teff (|Govev et al.l 120101 iRoias-Avala et all 120101 120121 
T12,M13a, M13b, N14, M14). Calibrations of these in¬ 
dices provide an efficient and simple method for mea¬ 
suring spectral type and temperature for the M dwarfs 
in our catalog. The relatively large widths of the wave¬ 
length ranges used in these measurements makes them 
insensitive to small errors in the RV measurement. Eor 
the M dwarfs in our catalog, we measur ed the Tpfp- 
sensi t ive indices used in M13b and N14 (iCovev et al.l 
120101 IRoias-Avala et afl 120121 iMann et al.l l2013al) fsee 
Tabled]). 

To measure spectral type for the M dwarfs in our 
catalog, we applied the NIR calibration developed in 
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Figure 1. The six separate telluric-corrected orders of the IRTF-SpeX SXD spectrum of GJ 1103A, which is a typical M4.5V star from 
our survey. Regions with significant telluric contamination are shaded in gray. 


Table 1 

Parameters of Teff-Sensitive Indices^ 


Name 


Range 1 (A) Range 2 (A) Range 3 (A) 


- 10172.5 10967.5 - 11052.5 11587.5 - 11672.5 

- 14655.0 16615.0 - 16645.0 17945.0 - 17975.0 

- 20277.5 22747.5 - 22782.5 23367.5 - 23402.5 

- 12300.0 13130.0 - 13330.0 13310.0 - 13510.0 

- 16150.0 16800.0 - 17000.0 17600.0 - 17800.0 

- 20900.0 22350.0 - 22550.0 23600.0 - 23800.0 


f20]^y J Index 10087.5 

H Index 14625.0 

Mann et al. K Index 20242.5 

Mann et al. 12014) H9O-J 12100.0 

^ove^e^a^ H20-i7 15950.0 

Roias-Avala et al. (2012'l H 2 O-A' 20700.0 

" Value of index = ((Fi> - {F 2 »/((F 2 > - (Fg)) 

N14, which is based on t he H 2 O-K 2 index developed in 
I Roias-Avala et all (|2012ll . This calibration yields “NIR 
M sub-types” on a uniform scale for our stars. The pre¬ 
cise relationship between this NIR M sub-type and the 
establ ished Palomar/Michigan State Unive rsity (PMSU) 
sc ale (iR.eid et al.lll995HHawlev et anilQQfiH is considered 
in iRoias-Avala et al.l (|2012ll and N14. and is known to be 
sensitive to [Fe/H] for at least the early M dwarfs. We 
also include in our catalog the measured PMSU spectral 
types from the N14 PMSU spectral type calibration over 
the range in which it is valid (M1-M4). 

To check the consistency of our NIR spectral type mea¬ 
surements, we considered the overlap between our sample 
and that of N14, which includes 152 M dwarfs. We found 
no significant offset between the respective measures of 
spectral type, and a scatter of half a subtype, consistent 
with the scatter expected from N14. 

This NIR spectral type measurement was necessarily 
the first in our chain of measurements, as it defined the 
template for our measurement of RV. The Teff-sensitive 
features are sufficiently broad that they are insensitive 
to RV shifts in our R ^ 2000 spectra. Nonetheless, after 
measuring the RV we re-measured the spectral type in 
the rest frame of the star, and this is the value we report 
in our catalog. 

3.2. Radial Velocity 


We measured the RV of each order of each spectrum 
by fitting a Gaussian to the cross-correlation (in log-A) 
of the observed spectrum with a spectral template of 
similar spectral type. We used thr ee spectral templates 
from the IRTF Co ol Stars Library (jCushing et al.ll2005l : 
iRavner et al.ll200^ : for stars M4 and earlier we used the 
M1.5V star HD 36395, for stars M5 to M6 we used the 
M5V star G1 51, and for stars later than M6 we used the 
MOV star LHS 2065. For each observation, we report the 
weighted (by SNR) average RV from all orders, and this 
is value we used to shift all spectra to a common frame 
for further analysis. 

To verify our RV measurements, we compared our 
(baryce ntric-corrected) RV m easurements to those mea¬ 
sured in iChubak et al.l (I2012D and N14 . For the 33 stars 
in common with |Chubak et al.l 1201211 , we find a me¬ 
dian offset of -1-7.6 km s“^ with a standard deviation of 
8.7 km s“^. For 152 stars in common with N14, we find a 
median offset of -1-7.9 km s“^ with a standard deviation 
of 13.2 km s“^. Figure [2] shows the differences in RV 
between this work and the literature sources. Among 
our common targets with the N14 study, there were 
seven with differences greater than 40 km s“^. Two of 
these (2MASS J07362513-b0704431/GJ3454 and 2MASS 
J17074083-I-0722066/GJ 1210) are known visual bina¬ 
ries with sub-arcse cond separation (jPravdo et al.l 1201)5 
l.lanson et al.l 120141 1 . Given the consistency of the RV 
measurements of the majority of the common targets. 
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$ Chubak et al. (2012) 
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-100 -50 0 50 100 

R^Lit. (km s-i) 

Figure 2. Differenc es between the RVs measured in this work 
and those measured in IChubak et a nonii) and N14. The median 
differences are shown as dashed lines, and are both approximately 
8 km s-l. 

we suspect that these highly discrepant targets are ei¬ 
ther RV-variable or binary stars. We also measure a large 
RV for one star (2MASS J06052936-b6049231/LHS 1817, 
— 113.9 km s“^) with similar magnitud e but the opposite 
sign as N14 and [Shkolnik et al.l p012l l. 

Tha t the systematic offset from N14 and iChubak et akl 
(|2012f) would be similar is unsurprising, as N14 imple- 
m ent a ^ 2.6 km s~^ c orrection to match the RV system 
of iChubak et al.l (l2012fl . The origin of the ^7 — 8 km s“^ 
offset for our data is unclear; since we rely solely on the 
arc lamp wavelength calibrations, we are subject to dif¬ 
ferent and possibly worse instrumental systematics than 
N14, who use telluric spectra to refine their wavelength 
calibrations. A systematic issue with the wavelength cal¬ 
ibration of approximately 8 km s“^ is also supported by 
the magnitudes of the inter-ord er an d RV template dif¬ 
ferences, as discussed in Section [3^ 

I n order to obtain ba rycentric RVs in the same system 
as (jChubak et al.l[2(TT^ and N14, the median offset above 
(^7.6 km s“^) should be added to the RVs in our catalog. 

3.3. Effective Temperature and Related Parameters 

Due to the difficulty of modeling M dwarf atmospheres 
and spectra, mo del-dependent measure ments of Taff can 
syste matically differ by 100 K or more (jRainurohit et al.l 
120131 and references therein). Interferometric Tes mea¬ 
surements have therefore been crucial in defining the M 
dwarf Teff scale, and in parsing the rel ations between 
the fu ndamental parameters of M dwarfs (jBovaiian et al.l 
[Ml and references therein). However, there are few 
targets with spectral type M5 or later that are bright 
enough or have large enough angular diameters to sup¬ 
port precise interferometric radius measurements. To en¬ 
able stellar radius measurements without interferometry, 
multiple groups have constructed empirica l calibrations 
based on photometry (jBovaiian et al.ll2012tl , spectral in- 
dices (Ml3b). or ind ividual absorption feature strengths 
(| Newt on et ai.l[2M^ . 

We measured Tes for the M dwarfs in our catalog using 
the H 2 O indices as discussed in Section [3Tl We used the 
J, i7. Kg Teff calibrations developed in M13b over their 


well-calibrated range (e.g. 1.0 < H 2 O—AT < 1.37), which 
corresponds to 3300 < Tea < 4800 K. M13b provide em¬ 
pirical relations for stellar radius, luminosity, and mass 
as a function of Tea, calibrated on a set of stars with 
interferometrically measured parameters. We favor the 
use of the TCs-band index and calibrations due to the 
robustness of this index to small changes in the telluric 
correction (as discussed in Section [TT|) . 

We also measured Tea, stellar radius, and stellar lumi¬ 
nosity using the i7-band at omic feature strengths stud¬ 
ied in iNewton et al.l (|2015ll . We measure the pseudo¬ 
equivalent widths (EWs) of these features using the con¬ 
tinuum and feature definitions in their Table 1. They 
construct relations using the strengths and ratios of Mg 
and A1 lines in M dwarfs, empirically calibrated on stars 
with interferometrically measured parameters. These re¬ 
lations use a similar calibration set to that of Ml3b, so 
the range of validity is similar: 3200 < Teff < 4800 K. 

3.4. Composition 

We measured [Fe/H] and [M/H] using the empirical 
spectroscopic calibrations from T12, M13a, M14, and 
N14. These calibrations were all developed using data 
from the same spectrograph (IRTF-SpeX) in the same 
settings (SXD) as ours, so their recipes can be applied 
directly to our data. M13a present separate calibrations 
for features in the J, H, Kg bands that yield both [Fe/H] 
and [M/H] for early-mid M dwarfs (M1-M5). They de¬ 
rived these calibrations from a blind search of the most 
[Fe/H] and [M/H]-sensitive regions of the spectra in 112 
stars with known [Fe/H] and [M/H], and are based on the 
EWs of a small number (four or fewer) of spectral regions 
in each case. The calibrations are valid for K5-M5 dwarfs 
and over approximately —1.0 <[Ee/H] < -1-0.56. We re¬ 
port the results of these calibrations separately for each 
target in our catalog. 

For a similar range of targets in spectral type and 
[Fe/H], we also report the ATs-band [Fe/H] calibration 
developed in N14, which depends only on the EW of the 
TCg-band Na i doublet. N14 calibrated this relation us¬ 
ing 36 M dwarfs with known [Fe/H] , and it is similar 
in for m to the relations developed in (|Roias-Avala et aH 
l2010f) and (|Roias-Avala et ^l2012f) . It is valid for MI¬ 
MS dwarfs and for a range of —1.00 <[Fe/H]< -1-0.35. 

For the late M dwarfs (M5 or later), the only applica¬ 
ble published [Fe/H] calibration is that of (M14). The 
M14 calibration is based on two strong features in the 
ATs-band, the Ca i triplet and the N a i doublet, which 
were also u sed in the calibrat i on of iRoias-Avala et al.l 
(|2010f) and IRoias-Avala et ahl (|2012ll for early-type M 
dwarfs. This calibration is based on 44 M dwarfs with 
known metallicities. It is an extension of the binary cal¬ 
ibration technique described above, which extends the 
calibration set by including some early-|-late M dwarf 
pairs where the earlier-type M dwarf [Fe/H] has itself 
been measured using empirical calibrations. This late- 
type calibration is valid for M4.5-M9.5 dwarfs and for a 
range of —0.58 <[Fe/H]< -1-0.56. The summary of the 
parameters of this [Fe/H] calibration and those above is 
shown in Table [21 

For general comparison with other catalogs, and for 
cases below in which it is useful to have a single [Fe/H] 
measurement for a given star, we elect to use the ATg-band 
calibrations of M13a for M1-M5 (see Section lOl) . along 
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with the output of M14 for M5 and later. We also report 
the [M/H] values where available. For the purposes of 
presenting and vetting our catalog, we primarily consider 
[Fe/H] as this is the specific quantity addressed by the 
bulk of work on the topic of M dwarf compositions. 


3.5. Activity and gravity-sensitive lines 

Our spectra extend to approximately 0.8 /rm on the 
blue end, and so contain the 8200ANa i doublet (here¬ 
after Na i) and the 8600ACa ii triplet (hereafter Ca ii), 
strong features that have shown promise as indicators 
of ac tivity (e.g. iKafka fc Honevciittl 120061: [Barnes et al.l 
1201411 and surface gravity ( Schlieder et al.l 1201211 . Our 
measurements of these features for a subset of the ob¬ 
served M dwarfs showed a strong relationship with the 
Teff and [Fe/H] measures in these stars, and can be 
used to derive information about their radii, luminosities, 
and possib ly g-enrichment. T hese measurements are dis¬ 
cussed in ([Terrien et al.ll201^ : we repeat the salient de¬ 
tails here. 

We measured the EWs of the Na i and Ca ii features 
(EWNai , EWcaii ) for each of our targets. Eor Na i, we 
chose li ne and pseudo - conti n uum regions sim i lar to those 
used in iMartiii et al.l (1201011 , iSchlieder et al.l (|2012H , and 
iGalvez-Ortiz et al.l (l2014lb For Ca ii, we defined feature 
regions for each of the three lines, as well as pseudo¬ 
continuum regions blueward and redward of each line. In 
each case, we fit a line to the pseudo-continuum regions in 
order to define the pseudo-continuum, and calculated the 
feature EW relative to this line. Eor EWNai and EWcaii 
we follow the convention that positive EW corresponds 
to line absorption. 

We caution that the Ca ii and Na i features are em¬ 
bedded in regions of strong TiO and oth er absorption (as 
can be seen in e.g. iRavner et al.l ([200911 1. so the EWs we 
measure are likely strongly sensitive to effects of the line 
species themselves as well as to the species that define the 
local pseudo-continuum. Moreover the SNR in the rele¬ 
vant orders of the SpeX spectra (orders seven and eight) 
is typically only a third of that in the J, H, Kg bands. 
Accordingly, these feature strengths should be carefully 
vetted before they are used in any application. The re¬ 
gions we defined for this analysis are listed in Table [3] 
and are shown on example spectra in Figure jS] Figure 
[3] also shows the three (possibly young) M dwarfs in our 
sample for which the Ca ii lines are in emission. 


3.6. Estimation of Uncertainties 

In order to estimate our measurement uncertainties 
(for all quantities other than RV), we rely on a set of 21 
M dwarfs with multiple observations. We calculate the 
RMS of the epoch-to-epoch differences in each measure¬ 
ment to obtain an observational estimate of our internal 
uncertainty. These values are shown for [Fe/H], [M/H], 
M spectral sub-type, Na i, and Ca ii in Figure (H for 
the Tfiff -sensitive indices in 5l and for the i/-band fea¬ 
tures of iNewton et al.l ([2015 1 in Figure | 6 | For the Ca ii 
and Na i features, we exclude three pairs with very low 
SNR (< 60). The observed RMS differences naturally in¬ 
clude effects that arise from variable telluric conditions 
and by-hand choices in the telluric correction pipeline. 
The resultant variable quality of the telluric correction 
can induce correlated noise which dominates our errors 
at the high SNRs we achieve (N14). 


We note that the Teff-sensitive indices may be affected 
by the broad characteristics of the telluric correction, in¬ 
cluding the flux calibration (which uses SIMBAD B, V 
magnitudes. Section [2^ and the construction of the tel¬ 
luric spectrum. By repeating the telluric correction with 
varying B, V magnitudes for a small set of stars, we found 
that the J-band indices were generally much more sen¬ 
sitive to differences in B, V magnitude and small differ¬ 
ences in the by-hand components of the telluric correc¬ 
tion. We favor the use of the A's-band indices (H 2 O-K 2 
or the M13b ATs-index), which are stable to a few percent 
or better with B, V variations of 0.5 mag. 

The Teff-sensitive indices can be used to estimate stel¬ 
lar Teff (also mass and bolometric luminosity) as de¬ 
scribed in M13b, who present uncertainties of 70— 100 K 
for these relations. Our adopted measurement uncertain¬ 
ties for their Teff-sensitive indices ('^ 0.01 — 0.02) corre¬ 
spond to median uncertainties of 110 K (J), 170 K (H), 
and 78 K (Kg) when propagated through the Teff-index 
relations from Ml3b. 

The iJ-ban d atom ic feature strengths presented in 
INewton et al.l ([2015H can similarly be used in empiri¬ 
cal calibrations to derive stellar Teff (residual scatter 
= 73 K), luminosity (0.049 L©), and radius (0.027i?©). 
Our repeated measurements of these i7-band features 
(Figured]) suggests errors of 0.1 — 0.2 A, which corre¬ 
sponds to measurement errors of approximately 90 K, 
0.04i? 0 , and 0.07 logL©. 

To estimate the RV measurement uncertainty, we re¬ 
port the standard deviation of the RVs measured from 
each spectral order. This empirical measure of uncer¬ 
tainty is conveniently available for all targets, and avoids 
potential issues of RV variability which may occur us¬ 
ing the epoch-to-epoch RMS differences. The median 
RV spread is ~ 5.7 km s“^. We note that a significant 
component of this RV spread is contributed by system¬ 
atic offsets of ^ 2 — 15 km s—1 between spectral orders. 
We found these shifts to be variable depending on the 
choice of RV template spectrum, suggesting issues with 
the wavelength calibration, line shapes, or PSF changes. 
We also found overall RV differences of a similar magni¬ 
tude (±20 km s“^) for each star depending on the tem¬ 
plate used, suggesting that our RV precision is dominated 
by these systematic issues. This also suggests that the 
^ 7 km s“^ offset between our RVs and those in the 
literature could be an artifact of the same systematic 
wavelength calibration issue. 

Alternatively, it is possible to calculate measurement 
uncertainties based on the errors reported by the SpeX- 
Tool pipeline, which reports a combined noise that in¬ 
cludes photon noise, read noise, and residual background 
noise components. Measurement uncertainties can then 
be estimated using Monte Carlo realizations of the noise. 
This process does not naturally include correlated noise 
effects like those arising from the telluric correction, al¬ 
though an ad-hoc technique is presented in N14 that 
mimics the effects of correlated noise. Briefly, this pro¬ 
cedure involves convolving the noise “spectrum” with a 
1.5 pixel gaussian, which is a good approximation for 
the actual autocorrelation function for observed spectra. 
This procedure effectively spreads the noise for each pixel 
over the neighboring pixels, mimicking the effect of the 
systematic correlated errors. Applying this technique to 
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Table 2 

Parameters of primary empirical calibrations 


Calibration 

SpT Range 

Metallicity Range 

Error 

Mann et al. 

Mann et al. 

Mann et al. 

Newton et al 

Newton et al 

2013a) [Fe/H]^ 
'2013a) [M/H]k 
' 2014) [Fe/H]x 
~7M4'| [Fe/H] 

. NIR SpT 

K5.0-M5.0 

K5.0-M5.0 

M4.5-M9.5 

M1.0-M5.0 

M1.0-M9.0 

-1.04 < [Fe/H 
-0.70 < [M/H 
-0.58 < [Fe/H 
-1.00 < [Fe/H 

< -1-0.56 

< +0.50 

< +0.56 

< +0.35 

0.11 dex 
0.10 dex 
0.07 dex 
0.12 dex 

0.5 subtype 



21410461+1935353 
02175601+1225266 
04472104+2808531 
GJ 699 


21410461+1935353 
02175601+1225266 
04472104+2808531 
GJ 699 


21410461+1935353 
02175601+1225266 
04472104+2808531 
GJ 699 


2.10 2.15 2.20 2.25 2.30 


A(/rrn) 


Figure 3. Three stars in the survey sample that have strong emission lines, along with GJ 699 as an example of a typical M dwarf. Three 
spectral regions with strong emission lines are shown; the top highlights the NIR Ca triplet, the middle highlights the H Pa /3 line, and 
the bottom highlights the H Br 7 line. The top panel also identifies the Na I and Ca II features described in Section 13.51 


Table 3 

Definitions for 8200A Na I and 86 OOA Ca II 


Feature 

Center 

(A) 

Width 

(A) 

PC Blue^ 

(A) 

PC Red'^ 

(A) 

Nal 

8192 

34.0 

8149 - 8169 

8230 - 8270 

Ca III 

8500 

12 

8456 - 8488 

8512 - 8536 

Ca II 2 

8544 

12 

8512 - 8536 

8553 - 8573 

Ca II 3 

8664 

12 

8635 - 8652 

8672 - 8688 


PC = Pseudo-continuum region definition. 


our measurements yields error estimates that are con¬ 
sistently smaller than the observed epoch-to-epoch RMS 
differences, so we elect to report the more conservative 
RMS numbers as our estimated measurement uncertain¬ 
ties. 

We note that the EWNai and EWcaii are 
thou ght to be sensitiv e to c hromospheric activ¬ 
ity (iKa fka fc H onevcuttJ 120061 : iBarnes et all 12014 
iTerrien et al.l 120151 1. and may be variable on the 
timescales (days-months) of our repeated observations. 
The epoch-to-epoch RMS differences may therefore be 
an overestimate of our measurement error for these 
features. We also estimate a lower limit of 0.1 A on 
the measurement uncertainties for these features using 


a 100-trial Monte Carlo technique based on the errors 
reported by the SpeXTool pipeline. Following N14, we 
include an ad hoc correction to account for the effects 
of correlated noise. The true measurement errors for 
Na I and Ca ii likely lie in the range bracketed by these 
estimates, 0.1 — 0.5 A. 

4. CATALOG DESCRIPTION 

We provide here a brief description of the spectra of 
the M dwarfs observed in this survey, as well as the set 
of measured and compiled values for each M dwarf. 

4.1. The Spectra 

The spectra for each of the 886 M dwarfs are available 
as a gzipped archive of standard FITS tables hosted 
at the Penn State Scholarsphere data hosting serviaHI. 
The file for each spectrum is identified by the 2MASS 
ID of the target star. The FITS headers contain the 
standard information output by the SpeXTool pipeline, 
including observation meta-data (instrument configura¬ 
tion, airmass, observing time, etc.) and the SpeXTool- 
generated file history. We append to this the measured 

^ https://scholarsphere.psu.edu/files/5712mq85s#.VWNMS9rBzGc 
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EW«.,i (A) EWc..i (A) 



Figure 4. For targets we observed multiple times, the comparison of the [Fe/H] (including our primary [Fe/H] as discussed in Section l3.4l 
and using the technique of N14), [M/H], M subtype, and Ca II and Na I feature strength measurements, from two different epochs. The 
black lines indicate a one-to-one relation. We adopt the observed RMS scatter in each case as our measurement uncertainty. 
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Figure 5. For targets we observed multiple times, the comparison of the Teff-sensitive index measurements from two different epochs. 
N14 refers to the bands defined in N14, and M13 refers to those defined in M13b. The black lines indicate a one-to-one relation. We adopt 
the observed RMS scatter in each case as our measurement uncertainty. 
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Figure 6. For targets we observed multiple times, the comparison of the lNewton et ahl l|2Q15l') feature strength measurements from two 
different epochs. The black lines indicate a one-to-one relation. We adopt the observed RMS scatter in each case as our measurement 
uncertainty. 
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RV shifts and barycentric correction for each spectrum, 
which can be used to velocity shift the spectra into a 
common reference frame. 

The spectra themselves are in the standard for- 
mat described fully in the SpeXTool documentation 
(jCushing et al.l[^004r FI. Briefly, each spectrum is a three- 
dimensional array, with the following axes: 

1. Order number (3 — 8). 

2. Data type (wavelength, flux, or variance). Wave¬ 
length is in vacuum, and is in the observer’s refer¬ 
ence frame. 

3. Pixel number. 

4.2. The Catalog 

Table 2] describes the columns included in our catalog. 
There is one line per star, including astrometry and pho¬ 
tometry compiled from the literature, measurements of 
spectral features, and derived stellar parameters. The 
full machine-readable table is available in the online edi¬ 
tion of The Astrophysical Journal Supplement. 

5. GENERAL DISCUSSION 

The overall distribution of our sample in NIR spectral 
type and [Fe/H] is shown in Figure [T] The [Fe/H] dis¬ 
tribution of our sample is centered around solar, with a 
median [Fe/H] of -1-0.05 dex and a scatter of 0.23 dex. 
This sample is neither volume nor magnitude-limited, 
but is consistent with the quoted [Fe/H] measurement 
uncertainty (0.07 — 0.11 de x) and volume-limited mea- 
surements of the FGK stars ( Casagrande et al.ll2011ll and 
M dwarfs ([Johnson fc AddsII 2009I1 in the solar neighbor¬ 
hood, which find a median [Fe/H] of ^ —0.05 dex and 
a scatter of ~ 0.2 dex. The distribution of NIR spectral 
types, which peaks around M4 and has a long tail toward 
earlier spectral types, is mostly a byproduct of our ob¬ 
servational strategy, although the mass function of the 
Solar neighborhood also peaks around this spectral type 
([Henry et al.ll2006ll . 

5.1. Internal Comparisons 

The majority of M dwarfs in our catalog have spectral 
types of M5 or earlier, and so can have their [Fe/H] mea¬ 
sured with multiple techniques. We focus on the most 
recent of these calibrations (from M13a and N14), along 
with the iJ-band calibration of T12 for comparison. Our 
measures of [Fe/H] are compared to each other in Fig¬ 
ure m which clearly shows systematic offsets among the 
various calibrations of up to 0.2 — 0.3 dex. In general, 
the H and J-band techniques are in agreement, while 
the A's-band techniques show good agreement up to the 
saturation of the N14 calibration around 0.3 dex. 

To explore whether issues with the wavelength calibra¬ 
tion could be responsible for these systematic differences, 
we tested the sensitivity of the various [Fe/H] calibrations 
to RV shifts. We measured [Fe/H] for each target 50 
times with a normal distribution of RVs centered around 
the original measured RV and with a spread equivalent 
to our estimated systematic and random measurement 

® Up-to-date manuals are available at the IRTF-SpeX website: 
http://irtfweb.ifa.hawaii.edu/~ spex/ 



NIR M Spectral Type Normalized 

Density 


Figure 7. The spectral type and [Fe/H] distributions for all M 
dwarfs observed in our survey. Our [Fe/H] values are the average of 
the J, H, iCs-band calibrations from M13a for targets with spectral 
type M5 and earlier, and the [Fe/H] derived from Ks band features 
using the technique from M14 for targets with spectral type later 
than M5. The hatched region shows [Fe/H] values that exceed 
the calibrated range. The projection plots show a kernel density 
estimation for each parameter, clarifying the peaks, widths, and 
tails of each distribution. 

errors (typically 10-20 km s“^). Among the M13a cali¬ 
brations, we found that the mean scatter in the Ag-band 
calibration was best (at ~ 0.02 dex), compared to the 
A-band (^^ 0.03 dex) and J-band (^ 0.04 dex). This 
suggests that the features chosen for the M13a J and 
A-band calibrations are marginally more sensitive than 
the Ag-band features to small wavelength shifts, and the 
magnitude of the expected RV errors is not sufficient to 
explain the observed offsets. We found that an RV shift 
of > 100 km s“^ can remove the offset between the Ag 
and A-band measurements, although this causes the J- 
band measurements to diverge further. In any case, such 
a large shift is not justified given the agreement among 
our per-order RV measurements, those for different RV 
templates, and between our overall RV measurements 
and literature values. 

The differences among our implementations of the 
M13a techniques may be due to undiagnosed issues with 
the wavelength calibration or other unidentified differ¬ 
ences between our reduction procedure and those of 
M13a (e.g. the telluric correction). We present measure¬ 
ments from each band in our catalog, and justify below 
the choice of the Ag-band calibration of M13a as our pre¬ 
ferred abundance calibration for M1-M5 dwarfs (Section 

ESI). 

5.2. Comparison to Other Catalogs 

To check our results and to probe the relationships 
between different measures of metallicity, we compared 
our values (using the calibrations of M13a, M14, and 
N14) to those presented by a number of re cent works 
([Roias-Avala et apl2(I12l: iGaidos et al.l 12014 N14) for 
common targets. iRoias-Avala et al.l ([20lll rely on the 
strengths of the Ag-band Na i and Ca i features along 
with the H 2 O-K 2 index to account of Tea; N14 rely 
on the strength of the Ag-band Na i feature alone; 
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Table 4 

M Dwarf Catalog Column Descriptions® 


Name 


Description 


Units 


ID 

RADEG 

DEDEG 

JMAG 

HMAG 

KSMAG 

E.JMAG 

E.HMAG 

E.KSMAG 

2MQF 

AMAG 

E.AMAG 

PMRA 

PMDEG 

R.PM 

PLX 

E.PLX 

R.PLX 

VMAG 

RV 

E.RV 

BC 

MH20J 

CH20H 

RH20K2 

MJI 

MHI 

MKI 

NSPT 

NSPT.PMSU 

MTEFFJ 

MTEFFH 

MTEFFK 

MFEHJ 

MFEHH 

MFEHK 

MMHJ 

MMHH 

MMHK 

MFEHL 

NFEH 

T12FEHH 

EWNA 

EWCA 

NMG148 

NMG150 

NK151 

NMG157 

NSI158 

NG0161 

NG0162 

NAL167A 

NAL167B 

NMG171 

N15TEFF 

N15RAD 

N15LOGL 

NOTES 


2MASS Identifier 

Right Ascension in decimal degrees (J2000) 

Declination in decimal degrees (J2000) 

2MASS J magnitude 

2MASS H magnitude 

2MASS Ks magnitude 

Error in 2MASS J magnitude 

Error in 2MASS H magnitude 

Error in 2MASS it’s magnitude 

2MASS quality flag 

UGAG4 A magnitude 

Error in UGAG4 A magnitude 

Proper motion in Right Ascension 

Proper motion in Declincation 

Reference for Proper Motion 

Parallax 

Error in parallax 
Reference for parallax 
V magnitude from UGAG4 
Radial Velocity 

Error in Radial Velocity (scatter among spectral orders) 

Barycentric Correction 
H 2 O-J index defined by Mann et al. 2013a 
H 2 O-H index defined by Covey et al. 2010 
H20-K2 index defined by Rojas-Ayala et al. 2012 
J Index defined by Mann et al. 2013b 
H index defined by Mann et al. 2013b 
K index defined by Mann et al. 2013b 
M Spectral Subtype defined in Newton et al. 2014 
PMSU Spectral Type defined in Newton et al. 2014 
Effective Temperature from J index defined in Mann et al. 2013 
Effective Temperature from H index defined in Mann et al. 2013 
Effective Temperature from K index defined in Mann et al. 2013 
Fe/H from J-band defined in Mann et al. 2013 

Fe/H from i7-band defined in Mann et al. 2013 

Fe/H from R-band defined in Mann et al. 2013 

M/H] from J-band defined in Mann et al. 2013 
M/H from iJ-band defined in Mann et al. 2013 

M/H from R-band defined in Mann et al. 2013 

Fe/H for late-type M dwarfs defined in Mann et al. 2014 
Fe/H defined in Newton et al. 2014 
Fe/H from iJ-band defined in Terrien et al. 2012 
Equivalent Width of 820nm Na I doublet 
Equivlanet Width of 860nm Ca II triplet 

Equivalent Width of 1.48 /rm Mg feature defined in Newton et al. 2015 
Equivalent Width of 1.50 /rm Mg feature defined in Newton et al. 2015 
Equivalent Width of 1.51 /rm K feature feature defined in Newton et al. 2015 
Equivalent Width of 1.57 /rm Mg feature feature defined in Newton et al. 2015 
Equivalent Width of 1.58 /rm Si feature feature defined in Newton et al. 2015 
Equivalent Width of 1.61 /rm CO feature feature defined in Newton et al. 2015 
Equivalent Width of 1.62 /rm CO feature feature defined in Newton et al. 2015 
Equivalent Width of 1.67 /rm Al-a feature feature defined in Newton et al. 2015 
Equivalent Width of 1.67 /rm Al-b feature feature defined in Newton et al. 2015 
Equivalent Width of 1.71 /rm Mg feature feature defined in Newton et al. 2015 
Effective Temperature from Newton et al. 2015 calibration 
Stellar Radius from Newton et al. 2015 calibration 
Logarithm of stellar luminosity from Newton et al. 2015 calibration 
Notes 


deg 

deg 

mag 

mag 

mag 

mag 

mag 

mag 

mag 

mag 

mas/yr 

mas/yr 

mas 

mas 

mag 

m/s 

m/s 

m/s 


K 

K 

K 


O.lnm 

O.lnm 

O.lnm 

O.lnm 

O.lnm 

O.lnm 

O.lnm 

O.lnm 

O.lnm 

O.lnm 

O.lnm 

O.lnm 

K 

solRad 

[solLum] 


Table 4 is published in its entirety in the electronic edition of The Astrophysical Journal Supplement 


iGaidos et al.l (|2014D used the F-band spectroscopic cali¬ 
bration presented by M13a. 

The comparison of our values to those from other 
works is shown in Figure IHl along with the median 
and scatter in each case. Our application of the M13a 
TAg-band [Fe/H] calibration shows good agreement with 
these other works, although we find systematically higher 
[Fe/H] for those stars with the highest [Fe/H] values com- 
pared to values reported in lB.oias-Avala etldl (|2012[ ) and 
(N14). This behavior may be a byproduct of the fea¬ 


tures used in each case, whose sensitivity to [Fe/H] can 
saturate (as noted in N1 4). Particularly, th e Na a nd Ca 
ATg-band features used in iR.oias-Avala et al.l (|2012f ). N14, 
and M14 are known to saturate as [Fe/H] approaches 0.5 
dex. 

The H and J-band calibrations are known to have 
larger scatter (M13a), and this is borne out in our com¬ 
parisons of these measurements to other works. In both 
cases, the [Fe/H] we measure for stars with Solar [Fe/H] 
and higher is systematically 0.1 —0.2 dex lower than that 
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Figure 8. Comparisons between the different measures of [Fe/H] included in our catalog (vertical axis, identified within each subplot), for 
all M dwarfs observed. The systematic offsets between the iiTg-band measurements and the J, H-band measurements are evident. Each plot 
lists the median offset (A) and standard deviation (a) of the differences, although these do not capture higher-order systematic differences. 


reported in other works. This offset, which is not present 
for the it's-band [Fe/H] values, may result from structure 
remaining in the residuals of these calibrations as shown 
in e.g. Figure 5 of M13a. 

There are several M dwarfs for which we measure 
later spectal type than M5, which also have [Fe/H] re- 
ported in other works. For these stars, the N14 and 
iRoias-Avala et all (|2012D calibrations appear to perform 
well compared to the M14 calibration, likely due to their 
common reliance on the i^s-band Na i feature. 

Our application of the N14 [Fe/H] calibration is con¬ 
sistent with the published [Fe/H] measurements in their 
work, ver ifying our implementation . For targets in com¬ 
mon with lRoias-Avala et al.l (|2012D . our [Fe/H] measure¬ 
ments are systematically higher, consistent with the off¬ 
set for later spectral types observed in Figure 18 of N14. 

A similar catalog to the one presented here is that be¬ 
ing dev eloped for the CAR MENES M dwarf exoplanet 
survey (iQuirrenbach et ^ 120141) , for which an exten- 
sive optical spectroscopic su rvey has been carried out 
(jAlonso-Floriano et al.l]^15[ ). There are 276 targets in 
the CARMENES catalog with coordinate matches in our 
catalog to within 5”. We plot the respective measures 


of composition and spect ral type in Figure HOI As an 
indication of abundance, lAlonso-Floriano et al.l 1)20151 ) 
measure the [Fe/H]-sensitive ( para meter ('based on TiO 
and CaH absorpt ion in the optical, iLeoine et all l2007t 
iDhital et al.ll20i^ . We find that C is weakly correlated 
with our measured [Fe/H] (Kendall’s r = 0.25, two-sided 
p = 1.3 X 10“®). This is unsurprising since is a coarse 
indicator of abundance, which has been primarily used 
for discriminating between dwarfs and subdwarfs. For 
all targets in the overlap sample , our [Fe/H] measure¬ 
ments and the C measurements of lAlonso-Floriano et all 
(|2015l) are both consis tent with dwarf compo s itions . To 
measure spectral type, lAlonso-Floriano et al.l (|2015D use 
a suite of spectral standard template matching and spec¬ 
tral index based techniques, and report spectral types 
to 0.5 subtypes. O ur spectral type estim a tes b roadly 
agree with those of lAlonso-Floriano et al.l (120151) : our 
NIR M subtype measurements have a median offset of 
0.37 with a residual standard deviation of 0.52, while 
our PMSU spectral types have a median offset of 0.12 
subtypes and a residual st andard deviation of 0.41. 
lAlonso-Floriano et al.l ()2015D use PMSU standards for 
much of their reference set, so it is expected that the 
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Figure 9. For stars i n common with o ther recent surveys 
IIRoias-Avala et al.l 12013 : IGaidos et al.( 120141 . N14), the compari¬ 
son to [Fe/H] values we measure. Each plot lists the median offset 
(A) and standard deviation (a) of the differences, although these 
do not capture higher-order systematic differences. 



Fig ure 10. For stars in comm on the CARMENES target cata¬ 
log (lAlonso-Floriano et al.|[2015ll . the respective measures spectral 
type. 

PMSU measurements should show better agreement. 


5.3. Preferred Measure of [Fe/H] for Early-Mid M 
Dwarfs 

In light of the systematic offsets shown in Sections 15.II 
and 15.21 it is useful to choose and justify a preferred 
measure of [Fe/H] for M1-M5 dwarfs. We choose the 
/Lg-band calibration of M13a for the following reasons: 

• The breadth and number of calibrators ex- 
ceeds those of other s imilar calibrations (T12, 
iRoias-Avala et al1l2012l : iNewton et al.ll201^ . 


The better agreement between the ATg-band mea¬ 
surements and other literature measurements of the 
star lSection l5.2|) . which are themselves b ased on 
ATg-band features (IRoias-A vala et a ni2012l , N14) 
and visible features (|Gaidos et al.11201^ . 


• The stability against small RV shifts of the ATg- 
band measurements compared to those of the J 
and H-band of M13a. 


5.4. [Fe/H] and Stellar Kinematics 

The stellar population of the Galaxy comprises multi¬ 
ple sub-populations with distinct kinematic, spatial, and 
compositional characteristics. These populations are of¬ 
ten grouped into three components—the thin disk, the 
thick disk, and the Galactic halo (although the disk pop¬ 
ulations may be better-represented by a continuous dis- 
tribution of popu lations than a bi-modal distribution, 
iBovv et ahl I20I2I1 . Large scale surveys of FGK stars 
and luminous evolved stars have found rich information 
about the history and evolution of the Galaxy within 
the detailed abundances of t hese different populations 
(e.g. iBensbv et al.ll20(3l. I2014lj. Similar large-scale stud- 
ies of M dwarfs ('e.g. iWest et al.l 120111 : iBochanski et al.1 
[Ml have also been carried out, and measurements of 
metallicity-sensitive band indices have enabled the iden¬ 
tification and study o f the oldest metal-poo r M subd¬ 
warf populations (e.g. IBochanski et al.ll2007f) . Refined 
metallicity measurement techniques are under develop¬ 
ment, and have demonstrated great potent ial for appli- 
cation with these large M dwarf catalogs (iDhital et al.1 
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[Mnl [Ml). Precise metallicity measurements of the 
M dwarfs in the solar neighborhood reveal correlations 
between composition and stellar kinematic and activity 
characteristics that are consistent with the increasing age 
and decrea sing level of activity from the thin disk to the 
thick disk (jRoias-Avala et aLfl2012tl . These results sug¬ 
gest the promise of precise abundance measurements in 
these long-lived stars. 

We explored whether the [Fe/H] values in were consis¬ 
tent with the expected trends, as a check on our [Fe/H] 
values and to probe for outliers and systematics. We cal¬ 
culated Galactocentric UVW motions for all of our tar¬ 
gets with parallax measurements. Combining these par¬ 
allaxes with proper motion values from LGll and LSPM- 
N and RV measurements from the literature where avail¬ 
able (or our RV measurements where literature RVs were 
unavailable), we have 283 targets with relatively precise 
UVW velocities (parallax error < 10%, literature RV or 
measured RV error < 10 km s“^, no indications of bina- 
rity). These velocities are shown in the Toomre diagram 
in the left panel of Figure 1111 along with ro ugh demar¬ 
catio ns of the thin disk, thick disk, and halo (jFuhrmannl 
1200411 . The Toomre diagram is effectively a projection 
of the rotational kinetic energy of each star relative to 
the local standard of rest, with stars of hi gher energies 
at increasing radii (jSandage fc Foutjll987[l . Consistent 
with expectations from galactic kinematics, those targets 
with the highest [Fe/H] values are kinematically within 
the thin disk region, and targets with lower [Fe/H] val¬ 
ues are spread across the entire range of velocities. The 
right panel of Figure |TT] shows the total Galactocentric 
velocity plotted against H a EW, a common activity in¬ 
dicator for M dwarfs (e.g. iReid et al.l 119951 : iWest et HI 
12004 IWalkowicz fc Hawlevll2009[) . The concentration of 
chromospherically active M dwarfs at low Galactocen¬ 
tric velocities, compared with the absence of active M 
dwarfs at higher velocities, suggests that we are indeed 
probing both the thick and thin disks, consistent with 
the chemical and kinematic properties of the M dwarfs 
in our catalog. The boundary be tween the thick an d thin 
disk populations is not discrete (llvezic et al.ll2008ll . and 
there are indications that the compositional gradient be- 
tween the thick and th in disks should be less than 1 dex 
(jBochanski et al.ll2007[l . consistent with Figure ITT] 

The most significant outlier in these plots is GM Dra 
AB, a well-st udied short-period ec li psing binary M dwarf 
system le.g. [Morales et 'ah! 120091: iTerrien et al.l \2012ak 
iFeiden fc Chabovei[l2014ll . Constraints fr om the cooling 
age of a nearby co-moving wh ite dwarf ([Morales et al.l 
l2009t IFeiden fc Chabov^ (201411 support the old age of 
this system, consistent with it s high velocity, low [Fe/H] , 
and possible g-en richment (IFeiden fc Chaboveii 12014 
ITerrien et al.l [201511 . That this system is chromospher¬ 
ically active at its old age is not surprising due to its 
short period (1.27 d) and the tidal interactions that have 
presumably taken place. 

5.5. Metallicity and Stellar Colors 

The infrared colors of M dwarfs have long been known 
to show sensitivity to stellar metallicity. This is in 
large part due to two primary sources of continuum 
opacity: H“ absorption and H 2 O absorption. Lower- 
metallicity stars have higher H~ and H 2 O absorption, 
leading to a loss of flux especially in the H-band 


([Mouldl Il976t iLeggetd 1199211 . The effects of metallic¬ 
ity on infrared colors have been exploited by many au- 
thors to construct p hotometric metallicity calibrations 
([Johnson et al.l 12014 M13a, N14), which are comple¬ 
mentary (observationally cheaper but less precise) to the 
spectroscopic [Fe/H] provided in this work. 

As another check of the integrity of our catalog of 
metallicities, we show in Figure [T^ the {J — H, H — Kg) 
colors along with the [Fe/H] values for each of our tar¬ 
gets. We confirmed that our data show a tight rela¬ 
tion between position in the {J — H,H — AT^j-space and 
[Fe/H], a connection established throu gh exaininatio n 
of kinematically-defined populations in iLeggettI (1199211 . 

Also shown ar e the kinematic age boundaries derived in 
iLeggetd ([199211 . transformed to the 2MASS photometric 
system using the updatec0 2MASS color transformations 
([CarDentei[l2001ll . 

5.6. Planet Hosts 

Stellar composition is of fundamental importance in 
studies of exoplanets. The properties of observational 
relationships relating t o composition, includin g the gi¬ 
ant planet-metallicity (iFischer fc Valenti] 1200511 and the 
planet mass-metallicitv (iGhezzi et al.ll2010ll correlations, 
provide unique insight into the environments and mech¬ 
anisms of planet formation. Moreover, reliable planet 
host star compositions are necessary in order to apply 
empirical relations hips to derive the stellar luminosities 
and temperatures (iBovaiian et al.l[2012l : iMuirhead et al.l 
I 2 OI 4 iMann et al.ll2015[l. and hence the location of the 
Habitable Zone ('e.g. IMuirhead et aI1l20124 M13b). 

This breadth of impact motivates the empirical and 
theoretical efforts to define an M dwarf metallicity scale. 
Multiple groups have already shown that the giant 
planet- metallicity correlation extends to the M d warf 
regime ([Johnson fc ApdsI 120091: [Neves et all l2013bll . A 
well-calibrated and precise metallicity scale is crucial for 
these and future studies of M dwarf planet populations. 

The metallicity values in this catalog are derived from 
the most precise and well-calibrated scales available for 
M dwarfs, and it is valuable to compare our results with 
those of previous studies for planet-hosting M dwarfs. 

Our catalog includes 16 M dwarfs that are known to 
be planet hosts, and in Table [5] we show our measured 
[Fe/H], [M/H], and H 2 O-K 2 based spectral types. Gen¬ 
erally, we measure metallicities consistent with those 
of previous NIR spectroscopic studies. In many cases 
our [Fe/H] measurements are substantially different from 
previous photometric estimates, which is expected as 
those calibrations are inherently less precise and use 
fewer calibrators than the relations we employed. 

6 . CONCLUSION 

We have presented NIR spectra and spectroscopic pa¬ 
rameters for 886 nearby M dwarfs. This catalog rep¬ 
resents the largest such compilation of M dwarf abun¬ 
dances, Teff-sensitive indices, and NIR spectra yet pub¬ 
lished, and includes: 

• [Fe/H] and [M/H] measurements (uncertainty ~ 

0.1 dex) for M1-M5 dwarfs using J, H, and Kg- 
band features, derived from multiple published cal¬ 
ibrations. 

® http: //www. astro. caltech. edu/-jmc/2mass/v3/traiisf ormations/ 
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Figure 11. (Left): The Toomre diagram for survey targets with measured parallax. The halo, thick disk, and thin disk regions are 
denoted. There are no metal-rich M dwarfs with thick-disk kinematics in our sample, consistent with expectations. A typical error bar 
is plotted in the lower left. (Right): The total Galactocentric velocity versus Ho EW (positive indicates emission and therefore more 
choromospheric activity). Indicators of both [Fe/H] and activity (Ho) suggest that our catalog probes both the thick and thin disks. 


Table 5 

Measured parameters for known planet hosts 


Name 

B05‘^ 

[Fe/H] 

JAOO'^ 

[Fe/H] 

W09^ 

[Fe/H] 

SL10'‘ 

[Fe/H] 

RAIO^ 

[Fe/H] 

RA12'‘ 

[Fe/H] 

N12^ 

[Fe/H] 

Tl2a 

[Fe/H] 

N14^ 

[Fe/H] 

other 

[Fe/H]'= 

This Work 
[M/H]*^ NIR SpT 

GJ 1214 


-1-0.03 


+0.28 

+0.39 

+0.20 

+0.03 


+0.05 


+0.40 

+0.27 

M4.5 

HIP 57050 

-0.02 

-1-0.32 

-1-0.07 


+0.12 

+0.04 


+0.05 



+0.04 

-0.03 

M4.0 

GJ 876 

-1-0.03 

-1-0.37 

-1-0.02 

+0.23 

+0.43 

+0.19 

+0.14 

+0.11 

+0.14 


+0.31 

+0.19 

M3.5 

GJ 317 







+0.22 

+0.31 

+0.22 


+0.43 

+0.30 

M3.5 

GJ 581 

-0.25 

-0.10 

-0.10 

-0.22 

-0.02 

-0.10 

-0.18 

-0.09 

-0.20 


-0.02 

-0.03 

M3.0 

GJ 179 


-1-0.30 

-1-0.02 

+0.20 


+0.23 

+0.14 

+0.18 

+0.12 


+0.25 

+0.12 

M3.0 

HIP 79431 

-1-0.16 

-1-0.52 


+0.35 

+0.60 

+0.46 


+0.50 



+0.78*' 

+0.55 

M2.5 

GJ 436 

-0.03 

-1-0.25 

-0.05 

+0.10 

+0.00 

+0.05 

+0.01 

+0.02 

-0.03 


+0.00 

-0.03 

M2.5 

GJ 849 

-1-0.14 

-1-0.58 

-1-0.20 

+0.41 

+0.49 

+0.31 

+0.24 

+0.22 

+0.22 


+0.50 

+0.33 

M2.5 

GJ 3470 










+0.20*= 

+0.27 

+0.15 

M2.0 

GJ 15A 










-0.36f 

-0.28 

-0.17 

M2.0 

GJ 649 

-0.18 

-1-0.04 

-0.13 

-0.03 

+0.14 

-0.04 


-0.05 



+0.04 

-0.01 

Ml.5 

GJ 433 







-0.13 


-0.17 


-0.03 

-0.06 

Ml.5 

Kepler-138 










-0.18S 

-0.21 

-0.16 

Ml.O 

WASP-80 










-0.14'= 

+0.13 

+0.07 

M0.5 

WASP-43 










+0.01= 

+0.40 

+0.28 

M0.5 


^ BOSiBonfils et al.l II20 05I1. JA09 I Johns on & Appsl 11200^ 1. W09 Woolf et al.l <200911. SLlOISchlau fman fc Laughlm| l|2010ll. RAlO lRoias- Ayala et all 
1120101 ). H,A12 ir{.oiS^^vala et all II2012I) . N12 INeveseran II 20 IS ). T12 lTCrnenet al.l Il2012bl) . N 14 INeveseral.l 112014] ) 

^ Uncertainties in [Fe/H] are approximately 0.11 dex 
^ Uncertainties in [M/H] are approximately 0.10 dex 

This star is well outs ide the calibrated range (—1.04 < [Fe/H] < -10.56) of [Fe/H], 

Demorv et ahll2013l) 


I Neves et al.ll2013al) 


M/H] value from iMuirhead et al.l Il2012al ) 
lichen et al.l[2()Tll 


• [Fe/H] with similar uncertainties derived from Kg- 
band features for M5-M8 dwarfs. 

• Measurements of J, H,Kg-hsaid H 2 0-based and 
other indices which are known to track Teff. 

• Measurements of strengths of the 820 nm Na i dou¬ 
blet and the 860nm Ca ii triplet, which are sensitive 
to activity and surface gravity. 

• Stellar photometry and astrometry from the liter¬ 
ature. 


• RVs (uncertainty ^ 10 km s ^). 

We performed cross-checks among the different [Fe/H] 
calibrations within our catalog and in the literature. We 
uncovered some systematic issues with the implementa¬ 
tions of the J and /[/-band calibrations, which led us to 
select the /Ts-band calibration of M13a as the preferred 
measure of [Fe/H] for M1-M5 dwarfs where multiple cal¬ 
ibrations exist. Using this [Fe/H], we confirmed the ex¬ 
pected trends with Galactocentric velocity activity levels, 
and stellar colors. We also provide updated abundances 
measurements for many of the known M dwarf exoplanet 
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